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ABSTRACT 

The  Green's  function  Is  evaluated  for  the  Lundquist 
equations  linearized  about  uniform  magnetic  field,  density 
and  zero  flow  velocity.   It  is  further  assumed  that  all 
quantities  vary  with  two  space  variables  and  time  only  and 
that  the  constant  magnetic  field  lies  in  the  plane  of  the 
two  space  variables. 
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NYO-2886 

On  the  Green's  Function  for  Two -Dimensional 
Magnetohydrodynamic  Waves.  I 


I.   Introduction. 

In  this  paper  we  shall  be  concerned  with  wave  motion 
in  a  plasma  where  the  plasma  is  considered  as  a  continuous 

medium.   A  particularly  simple  set  of  equations  describing 

2 
such  a  medium  is  the  Lundquist  equations^  which  represent 

an  isentropic  fluid  with  scalar  pressure  tensor  and  infinite 
conductivity,  and  in  which  displacement  currents  have  been 
ignored.   We  shall  examine  only  small  disturbances  Imposed 
on  a  plasma  initially  with  zero  flow  velocity  and  constant 
matter  density  and  magnetic  field;  that  is,  we  shall  employ 
the  Lundquist  equations  linearized  about  these  quantities. 
As  further  restrictions  waves  depending  Just  on  the  x  and 
y  coordinates  will  be  included  and  the  unperturbed  magnetic 
field  is  taken  in  the  direction  of  the  x-axis.   For  such  wave 
motion  we  shall  give  the  Green's  function,  or  fundamental 
solution,  which  represents  the  wave  emanating  from  a  point 
source.   The  Alfven  and  longitudinal  waves  are  automatically 
decoupled  and  many  simplifications  result.   That  the  magnetic 
field  is  divergenceless  introduces  complications.   This  con- 
dition is  in  fact  an  initial  condition  rather  than  a  dynamical 
equation,  and  the  techniques  used  are  not  naturally  adapted 
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for  the  imposition  of  initial  conditions.   With  an  artifice 
we  shall  be  able  to  achieve  the  result  of  imposing  the  condi- 
tion without  introducing  inconsistencies.   We  shall  treat  the 
Green's  function  as  a  distribution^  and  determine  it  by- 
inverting  its  Fourier  transform  obtained  from  the  partial 
differential  equations.   The  interpretation  as  a  distribution 
is  essential  as  we  shall  have  to  evaluate  many  singular 
integrals  and  the  Green's  function  itself  is  not  integrable 
considered  as  an  ordinary  function. 

The  author  wishes  to  thank  Dr.  Clifford  Gardner  for 
numerous  stimulating  discussions  and  Dr.  Harold  Grad  for  his 
suggestions,  many  of  which  have  been  incorporated  into  the 
paper. 
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II.   Preliminary  Material. 

In  rationalized  MKS  units  we  may  write  the  equations  for 
the  fundamental  solution,  or  Green's  function,  for  hydro- 
magnetic  waves  derived  from  the  linearized  Lundquist  equations 
as : 


(la)      if  +  Po^-  ^=  S  6(r)5(t) 


(lb)   Po  1^  +  a^Vp  -  1  (7XB)XB^  =  s^  6(r)6(t) 


(Ic)      If  -  VMuXB^)  =  ^  5(r)6(t), 


where  p^,  B^,  a^,  and  \x     are  the  constant  mass  density, 

magnetic  field,  ordinary  sound  speed,  and  magnetic  permeability 

of  the  undisturbed  fluid,  and  p,  B,  and  if  are  the  perturbed 

density,  magnetic  field  and  fluid  velocity.   ^,,  and  s'   are 

u       a 

I 
arbitrary  constant  vectors  and  S   Is  a  constant.   Although 

P 

we  refer  to  "the"  Green's  function,  clearly  we  mean  expres- 
sions for  p,  u,  and  B  satisfying  (1).   We  shall  give  such 
expressions,  which  have  the  form  of  sums  of  various  functions 
multiplied  by  S^  and  by  components  of  S  and  S„  for  the 

p  U         C 

case  under  consideration.   Temporarily  we  write  the  solution 
In  terms  of  scalars,  vectors,  and  dyadlcs  as: 


PP   P     PU 


P  =  G..  S^  +  G..-  S^  +  GpB-  Sb 
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^   =  Sp  2p  +  ^Bu-  Su  +  Sb-  ^ 

While  the  system  (1)  is  perfectly  satisfactory  mathemati- 
cally, there  is  an  additional  equation  that  should  be  included 
so  that  the  quantities  become  physically  meaningful : 

(Id)      V  '   B  =  0. 

However,  from  (Ic),  we  infer 

^=  6(t)V(6(r)^3)  ^   0 

A  consideration  of  the  use  of  Green's  functions  immediately 
suggests  methods  to  avoid  the  contradiction. 

If  instead  of  the  delta  functions  on  the  right  hand 
sides  of  (la),  (lb),  (Ic)  there  were  the  functions  p(r,t), 
u(r,t),  B(r,t),  then  we  might  write  the  solution  of  the 
system,  which  includes  the  initial  value  problem  a^  a  special 
case,  as: 

p(?,t)  =  J  (dr')(dt')|Gpp(?-?',t-t')p(r',t') 

+  Gpy(r-r',t-t').  u(r',t')  +  Gp3(?-?' ,t-t ' ) .B(r ' ,t ' )| 
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u(r,t)  =  J  (dr')(df  )[^^p(r-?'-,t-f  )p(r',t' ) 


+  S'yy(r-r',t-t')-u(r',f  )  +  ?^(r-r' ,t-t' )  .B(r' ,f  ) 


— ^/-^ 


B(r,t)  = 


r 


(dr-  )(dt'  )[^Bp(^-^'  ,t-f  )p(r'  ,f  ) 


+  ^By(r-r',t-t').u(r',t')  +  G^gCr-r ' ,t-t' ) 'BCr' 


,f  )] 


For  such  a  system  to  be  physically  realizable  B(r,t)  must 

satisfy  V*B(r,t)  =  0.   We  shall  separate  GgB  Into  two 

parts,  ^BB  =  ^;bb  "'■  ^BB  '   ^°  t\\dit   for  any  ^(r,t)  with  zero 
divergence 


r 


0  = 


(dr')?i|^(r-r',t-t')-  t(r',f) 


'BB 


0=  V 


(dr')&)(^_-.^t-f  )-t(r',t'). 


'BB 


Thus  while  both  O^^  and  ^^|^  occur  In  the  Green's  function 
of  the  original  system,  for  physically  acceptable  problems  we 
may  Ignore  Gil^.   Further,  although  (Id)  is  not  satisfied 
With  the  full  S'gg,  with  only  %^      (Id)  is  satisfied. 
Obviously  such  a  separation  is  not  unique.   In  order  to  under- 
stand in  what  sense  the  separation  chosen  is  a  very  natural 
one  and  to  verify  its  desired  properties  we  consider  an 
alternate  form  for  (Ic): 


8  - 


(IC)     ll  -  Vx(ux^^)  =  7K(#^,6(r))6(t) 


-1 


Now  automatically  hV  -B/St  =  0.  Again  if  we  solve  the 
inhomogeneous  problem  with  sources  ■p(r,t),  u(r,t)  and 
B(r,t)  =  ^xX(r,t),  and  with  an  obvious  extension  of  notation; 


§(r,t)  =  r(dr')(dt')JGBp(r-r',t-f)p(r',t') 


+  ^^^{r-r^  ,t-t '  )  -uCr  •  ,t  •  )+  GgA^^"^'  '^"^  '  ^  ''^^^'  '^'^ 


and  the  other  equations  are  unchanged.  After  an  integration 
by  parts  we  shall  see  that  the  last  term  in  the  equation  for 
B  becomes 


J 


(dr'  )(dt'  )  G^J^r-r'  ,t-t'  )'t(r'  ,t-  ), 


'BB 


and  all  the  desired  properties  follow.   We  need  not  separate 

G„   or  Gt,   into  two  parts  as  they  contribute  an  automati- 
ap  jdu 

cally  divergence-free  magnetic  field.   Since  in  the  equation 
for  u  only  VxB  occurs,  only  the  divergence-free  part  of 
B  is  present  and  we  need  not  make  any  separation  here.   We 
do  not  replace  (Ic)  by  (Ic")  entirely  primarily  as  a  matter 
of  completeness.   However,  we  shall  finally  give  explicit 
results  only  for  GAg   and  not  for  G^g. 

All  quantities  depend  only  on  x  and  y  and  the  spatial 
delta  functions  in  (1)  now  represent  5(x)6(y).   The  Fourier 
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and  Laplace  transform  of  p(r,t)  is  defined  by 


r 


p(r,t)  = 


(2Tr)- 


(dl^)(da3)p(lf,a))e^(^*^-^^)  , 


where  the  oo  integration  nans  from  -«  to  +«  on  a  line 
Just  above  the  real  axis  in  the  complex  to  plane,  and  the 
ic  integration,  ^  =(k^,  k^,  0),   (dl^)  =  (dk^)(dky),  runs 
over  all  real  vectors   (l^x'  ^^ '      ^^  define  in  a  similar 
manner  the  transforms  of  u(r,t)  and  B(r,t).   The  symbols 
p,  u,  B   represent  interchangeably  the  functions  and  trans- 
forms; context  will  readily  determine  which  is  involved.   With 
the  further  definitions  and  substitutions 


r 


Bq  =  |Bq|(1,0,0)  =  iB^lb 


b2 
A_  =  r-4—  ,  the  Alfven  speed 


o   lap^ 


S  =  p  S 


_ii 


(2)    i   S  =  p  S 
^  '  u   ^o  u 


_ii 


Sb  =  I^o'Sb 


B  is  replaced  by   iB  |B 


p  is  replaced  by  p  p 
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the  equations   for  the   transformed  variables  become 

(3a)  oop   -  I?-u  =  ISp 

(5b)  cou   -  a^  itp   +  A^(^x(Bxl^))    =  IS^ 

(5c)  oM  +  k  x(ux'^)   =  1^ 

(5c')  aB+lcx(uxb)   =  -(kxS^). 

In  the   case  under  consideration  the  Alfven  wave  does 
not   couple  to  the  rest  of  the   system,   thus 


(4a)  ccu^  +  A^k^B^  =  IS^^ 


(4b)  cnB^  +  ^x  "z  =  ^\' 


and  u   and  B   do  not  occur  elsewhere  In  the  system.   The 
z        z 

equations  (4)  represent  an  Alfven  wave,  and  a  trivial  calcu- 
lation gives : 

(5a)     u^(r,t)  =6%if(S^  -  A^Sb  )6(x-A„t) 

L   z      z 

^^^u  +  Vb  )^(^  -^  V)f 

z      z  J 

(5t)      B^(r.,t)  =  iM  ((A^Sg  -  S^  )6(x-A^t) 

o  (^     z    z 


z     z  ^ 
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The  remaining  coupled  equations  give  the  longitudinal  wave. 
With  the  definition 

D  =  0)"^  -  a>V(a2  +  a2)  +  a^  a2  k^k^ 


the  solutions  of  the  equations  (5a),  (3b),  (3c)  for  p,  u  , 
(6a)      p  =  ^L(cd2-  A^k2)Sp  +  a3^(k.S^)-  k^A^S^ 

(6b)      u^  =  §fco(a>2-  A2)s^  cuk^a^dc  S  -  k  S   ) 

V  X  y      X 

+  k^(c.2-  k^A^la^Sp  +  a^A^  Vy(kyS  -  k^Sg  )) 

X     y  V 

(6c)      Uy  =  ^[a.(co2  -  k2a2)s^^+  Vy"^o\  "^  V^^o^p 

^  (^'-  ^x^o)Ao(^\-  ^x^By)] 

(6d)      B^  =  §  [a)(co2-  a2k2)S3  +  a^^   a^s^  +  co^k^S^ 

V  X  y 

(6e)     By  =  ^L{</-   a^k^lSg  -  "k^kya^Sp  -  a>\s^ 
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In  the  expressions  (6d)  and  (6e)  for  B^  and  B   the 
terms  in  the  second  bracket  with  the  leading  factor  i^'^g 
are  the  "unphysical"  part  of  the  Green's  function  forming 
G^^'.   The  remainder  of  the  terms  depending  on  Sg  form  ^gg  • 
To  Justify  the  separation  and  the  general  statements  above  let 
us  consider  equation  (5c')  instead  of  (3c),  so  that  we  replace 
Sg  by  i(kx'S»).   First,  all  the  terms  from  G^^^  vanish 
identically,  so  that  G^g^   contributes  nothing  to  a  physical 
problem.   If  we  define  j(r,t)  as  the  Fourier  transform  of 
icD(a)  -  a  k  )/D  then  exhibiting  explicitly  the  dependence  on 
S.   only,  we  have 


^x  =  Vly  J(r,t)  -.  .. 

2 


^  =  -\%^'^^'^^  ■"  '" 


For  the  inhomogeneous  problem 


B^(r,t)  = 


(dr')dt'  °-  j(r-r',t-t')A^(r',t)  + 


and  by  an  integration  by  parts  (with  no  surface  terms  as  J 
has  compact  support  for  all  t) 

/^  V    r  -  X     /^^      X  ^A  (r',t) 
B^(r,t)  =  (dr')dt'  J(r-r',t-t')  — +  ... 

^       J  ^y' 


=  J(dr') 


dt'  J(r-r' ,t-t' )B^(r' ,t' )  +  ... 
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and  similarly  for  B  (r,t) 

V 


By(r,t)  =  J(dr')dt'  J(r-r',t-f )By(r',t'). 


However,  we  note  that  these  are  precisely  the  integrals  that 
occur  if  we  include  oig  only.   Similar  remarks  apply  to 
B„(r,t),  but  we  need  not  even  consider  this  component  as  B_ 
does  not  appear  in  7*B  since  no  quantities  depend  on  z. 


III.   General  Results  on  the  Inversion  of  the  Transforms. 


Initially  we  shall  consider  a  general  transform  defined 


by; 


(7) 


F(r,t)  = 


(27r)- 


r 


do) 


J 


J 


Q(oD,k) 


where  Q(a),k:)   is  a  homogeneous  polynomial  of  degree  m 
which  is  the  characteristic  polynomial  of  a  hyperbolic  partial 
differential  equation  and  a,  3,  y  are  positive  integers 
whose  sum  is  m-2.   All  the  terms  in  equations  (6)  are  of 
this  form  except  that  a+p-py  =  n^-l>  however  by  the  relation 


rdcor(dk)e^(^-^-^^^(-ia)t(a3,k)  +  ikn(a),k) ) 


^ 


da)r(dk)ei(^-^-^^)t(a),k))+V(Jda3|(dk)ei(^-^-^^)^ 
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we  may  express  all  the  transforms  in  (6)  as  derivatives  of 
the  type  considered  in  (7).  We  shall  use  polar  coordinates 
in  the  k  plane  slightly  different  from  the  ordinary  ones; 
the  "magnitude"  of  k,  k,  will  range  from  -»  to  +00,  and  the 
polar  angle  6     from  0  to  tt;  the  range  -oo<k<0,  0<e<T 
corresponding  to  0<k<eo,   7r<0<2ir  in  the  ordinary 
usage.   In  such  a  system  (dk)=dk^dk  =  |k|dk  de  =  k  sgn(k)dk  de. 
The  symbol  k  =  k(0)  denotes  the  unit  vector  in  the  0   direc- 
tion, with  the  direction  9=0  to  be  specified  later. 

The  zeros  of  the  denominator  of  the  integrand  in  (7), 
considered  as  a  function  of  o),  are 


ojj  =  kVj(0)   J  =  1,  2,  3,  ...  ,  m, 


The  v,(0)  are  just  the  speeds  of  propagation  of  the  various 
modes  in  the  9   direction.   With  the  symbol  H(t)   for  the 
Heaviside  function,  and  after  the  co  integration,  we  have 

00  T 

(8)       P(r,t)  =-^^  [  sgn(k)dk  r  de(k  l^'Ck  )P 


A 


ik(k-r-v  (e)t) 
^^e ^  (vj(e))'v 

f^         Qjvj(e),  k(e)) 


where  Q^^(Vj(e),  k)  =  |^  Q(a3,  k) 


a3=Vj(e) 
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The  k  Integration  is  meaningful  as  a  distribution.   With  the 
integral 


dk  sgn(k)e^^^  =  ^ 


and  on  identification  of  the  appropriate  partial  fraction 
decomposition,  we  have 


v 


(9)       F(r,t)  =  5M  r 

27r^t  J 


de  — 


«(¥'  ^(«)) 


For  future  use  we  note  that  the  zeros  of  the  denominator  of 
(9)  correspond  to  unit  vectors  k  such  that 


k«r 


(10)      ^  =  Vj(k)    J  =  1,  2,  5,  ...,  m. 


A 


In  (8)  we  have  tacitly  assumed  that  the  roots  of  Q(co,k)  are 
simple,  however  we  may  derive  equation  (9)  easily  even  if  the 
roots  are  multiple  or  of  multiplicity  changing  with  9     by 
the  appropriate  modification  of  (8). 

Since  the  denominator  in  (9)  has  zeros  not  integrable  in 
the  ordinary  sense,  we  must  give  some  interpretation  to  the 
integral.   One  might  suppose  that  the  appropriate  interpreta- 
tion is  the  usual  finite  part  integral  considered  by  various 
authors,  e.g.,  Lighthill.^  In  that  case  the  integrand  in  (9) 
would  be  extended  into  the  complex  9   plane  and  the  Integral 
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would  be  taken  as  the  average  of  the  two  integrals  on  contours 
just  above  and  Just  below  the  real  axis.   Such  a  prescription 
follows  directly  from  the  more  standard  formulas.   However, 
normally  we  expect  finite  part  integrals,  which  are  indeed 
always  finite,  to  represent  ordinary  functions.   Here  we 
expect  the  answer  to  be  a  distribution,  which  in  the  ordinary 
sense  need  not  be  finite.   It  is  beyond  the  scope  of  this 
paper  to  derive  rigorously  the  correct  interpretation  of  such 
integrals;  we  shall  merely  describe  and  discuss  the  results. 
The  primary  thing  to  note  is  that  the  finite  path  or  contour 
integral  interpretation  is  indeed  valid  for  strictly  hyper- 
bolic equations,  i.e.,  when  the  wavespeeds,  v.(9), 
J  =  1,  2,  ...,  m,  are  distinct  for  each  value  of  9,  provided 
none  of  the  roots  v.(0)   is  of  the  form  a-k,  with  a  a 

J 

constant  vector.   For  the  non-strlctly  hyperbolic  case,  again 
with  no  speeds  linear  in  k,  we  find  that  F(r,t)   is  given 
by  the  contour  integral  plus  some  functions  discontinuous 
across  certain  planes  associated  with  the  speeds  that  are 
multiple.   For  hydromagnetic  waves  the  roots  are  multiple  but 
ultimately  the  associated  complications  cancel. 

In  passing  we  give  some  ideas  as  to  how  tlie  results  are 
derived  rigorously.   Provided  the  multiplicity  of  the  roots  is 
not  too  great,  we  may  write 
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IT 


m  /f,.^  „  /o\^-^^l-3 


^^^■OhW 


(k.r-vj(0)tr-'^  log|i^.r-Vj(9)t| 


J^  Qjv.(e),  1?) 


7 
which  follows  from  plane  wave  techniques  or  merely  on 

Integration  of  (9)  formally  several  times  in  x,  y,  and 

t.   By  taking  F(r,t)  as  the  real  part  of  the  integral  we 

should  like  to  remove  the  absolute  value  signs  from  the 

logarithm,  displace  the  contour  of  integration  and  then 

differentiate,  deriving  the  contour  integral  representation. 

All  the  steps  may  be  justified  in  the  strictly  hyperbolic 

A 

case  with  wave  speeds  not  linear  in  k.   If  the  equations 
are  not  strictly  hyperbolic,  but  without  linear  wave  speeds, 
we  find  that  there  are  complications  with  branch  cuts  of  the 
logarithm,  and  it  is  just  these  jumps  that  give  rise  to  the 
discontinuous  functions. 
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IV.  Longitudinal  Waves. 

For  longitudinal  waves  the  polynomial  Q(co,k)  in  (7) 
is  just  o)^-  cD^(a^  +  A^)k^  +  a^  A^  k\^  so  that  the  wave 
speeds  are  given  by 


(11) 


,^,a)  =  -v,(a)=fej-  |((a^A^)^.  Aa^A^(5j^)'/y 


v,(e)  and  v,(e)  being  the  fast  wave  and  V2(0)  and  V2^(0) 
the  slow  wave.   For  k^=  0,  Vg  =  v^^  =  0.   Thus  the  wave  speed 
locus  has  a  double  point.   Further  if  a^  =  A^  then  for 
!<•  =  +  1,  V,  =  Vp   and  v^  =  v^   so  that  another  double  point 
may  occur.   As  described  in  detail  elsewhere  ,  the  double  point 
in  the  characteristic  equation  introduce  certain  planes  in 


x,y,t  space  where  there  is  irregular  behavior  of  the  Green's 
function  not  given  by  the  principle  value  interpretation  of  (9). 
For  the  double  root  at  k  =  0,   the  plane  is  y=0, 
(x/t)^  <  a^A^  /  a^+A^,  and  when  a^  =  A^,  so  that  k^  =  1  1 
is  another  double  rootj  there  are  planes  x=±a^t,   (y/x)  <  l/4. 
However  for  the  combinations  of  terms  that  appear  in  (6),  the 
irregularities  associated  with  these  planes  cancel  Identically. 
In  other  words,  while  the  planes  in  question  are  possible 
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singularity  surfaces  for  systems  of  equations  with  this  char- 
acteristic polynomial,  for  this  particular  system  they  have 
no  effect  and  we  may  Ignore  them  and  write 


(12) 


t.t)  ^mi 


P(r,t) 


4Tr^t 


de  —^ 2_ 


c+c 


Q(%^,  k) 


with  the  understanding  that  In  (11)  equality  holds  Ignoring 
terms  of  no  effect  for  the  system  of  equations  for  hydromagnetlc 
waves.  In  (12)  C  Is  a  contour  In  the  complex  9   plane  from 
0  to  IT  Just  above  the  real  axis  and  C  Is  a  contour  Just 
below  the  real  axis  from  0  to  tt. 

If  we  now  choose  the  0=0  direction  as  the  vector  B 
and  call  the  angle  between  r  and  B  0,     and  finally 


introduce  the  auxllllary  variable  z  =  e 


21(0-0) 


we  find 


(15) 


F(r,t)  = 


2t'^H 


TT  r  1 


^ 


zdz 


C(z) 


where 


(U) 


.2. .2. 


!(z)  -  z  +  4z-^(l A   t  +    M —  e  ^j+  z  (6 ^ t 


.  S.1.1  ^^yl.  '-^t^.'^  e-^)  .  X 
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and  where  I   and  I    are  contours  just  inside  and  Just 
outside  the  unit  circle  taken  in  a  positive  sense.   Since 
a  +  p  +  Y  =  2  the  niomerator  is  even  in  k  and  when  multiplied 
by  z  and  considered  as  a  function  of  z  it  is  a  second 
degree  polynomial  in  z.  Thus  the  value  of  the  integral  is 
27ri  times  the  residue  at  poles  inside  the  unit  circle  minus 
the  residue  at  poles  outside  the  unit  circle.   Some  of  the 
geometrical  results  derived  below  are  well  known  and  follow 

Q 

from  general  properties  of  partial  differential  equations.-^ 
We  rederive  them  as  a  matter  of  completeness  and  to  show  how 
they  reappear  when  employing  transform  techniques. 

If  C(z)  =  0,  then  also  C(l7z)  =  0,   so  that  1/z  is 
also  a  root  of  C(z)  provided  zz"  ^  1,   or  z  lies  off  the 
unit  circle.   For  various  values  of  the  parameters  there 
are  three  possible  situations:   four  roots  of  C(z)  on  the 
unit  circle;  two  on  the  unit  circle,  one  Inside  and  one  out- 
side; or  two  inside  and  two  outside.   With  this  information 

A 

we  shall  now  work  with  the  variable  k  instead  of  z,   the 

A 

relation  being  real  values  of  k  correspond  to  real  angles 

A 

6     and  z  on  the  unit  circle,  and  complex  k  to  all  other 
values  of  z. 

As  was  stated  (see  (10)),  the  zeros  of  C(z)  occur 
whenever 


(15)      {tc-t)  -^  =  p^,   j  =  1,  2,  3,  4 


Vj(k) 
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where 


(16)      p/p2  =  t/t    , 


i.e.,  p  is  the  Inverse  point  in  the  unit  circle  to  r/t. 
To  exhibit  the  conditions  geometrically  let  us  consider 
Figure  1  representing  the  p  plane.   The  positive  abscissa 
corresponds  to  the  positive  x  direction,  and  the  positive 
ordinate  to  the  positive  y  direction.   In  terms  of  polar 
coordinates  {v,-\l/)     in  the  p  plane  we  plot  the  curves 
(l/vJ-^),   if),      j  =  1,2,5,^.   Condition  (15)  implies  (OA)(OB)  = 
(OC)   where  AC  is  perpendicular  to  OB.   Hence  OAC  and 
OCB  are  similar  triangles  and  OCB  is  a  right  angle.   To 
determine  the  real  k  for  which  (15)  holds  we  find  the 
intersections  v;ith  the  curve  l/v^  of  the  line  Ji     through 
p  and  perpendicular  to  the  line  from  p  to  the  origin. 
Intersections  at  infinity  must  not  be  neglected,  thus  lines 
parallel  to  the  ordinate  axis  have  two  intersections  at  infinity 
with  the   (l/v)   curve.   Clearly  any  line  must  intersect  the 
slow  curve  at  least  twice,  and  correspondingly  C(z)  always 
has  at  least  two  roots  on  the  unit  circle. 

Except  in  the  neighborhood  of  those  values  of  the  coef- 
ficients of  C(z)   for  which  there  are  multiple  roots,  the 
roots  will  be  analytic  functions  of  the  coefficients  and  hence 
of  r/t.   Except  at  these  special  points  the  residues  of  the 
integral  (13)  will  also  be  analytic  functions  of  r/t.   Since 
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there  is  never  more  than  one  root  of  C(z)  Inside  the  unit 
circle  or  one  root  outside, it  suffices  to  find  the  values  of 
"p  in  Figure  1  which  correspond  to  multiple  intersections  of 
th  line  ^(p)  with  the  l/v  locus.   If  £     has  a  multiple 
intersection,  then  l     is  tangent  to  the  curve  or  the  point 
un  the  curve  is  a  multiple  point.   But  the  multiple  points 
are  associated  with  real  angles  6     and  hence  do  not  affect 
the  integral.   The  points  in  the  p  plane  at  which  there  is 
the  possibility  of  non-regular  behavior  of  the  integral  {1J>) 
form  the  set  of  points  each  being  the  closest  to  the  origin  on 
a  tangent  to  the  curve  l/v.   But  by  definition  this  curve  is 
just  the  pedal  curve  of  the  original  curve.   The  inverse  in 
the  unit  circle  of  the  pedal  curve,  a  curve  in  the  ordinary 
r/t  space,  is  well  known  to  be  a  plane  section  of  the  ray  cone, 
the  ray  cone  being  a  surface  carrying  discontinuities  of  the 
Green's  function.   We  note  that  for  this  characteristic  poly- 
nomial certain  other  planes  may  also  carry  discontinuities 
not  present  for  this  system  of  equations. 

The  ray  cone  is  shown  in  Figure  2  when  ^  ^^  A  .   In 
regions  A  and  C  there  are  four  real  Intersections  of  the 
line   -2(p)  and  the  l/v  locus  and  in  region  B  tv/o.   In 
regions  A  and  C  there  are  four  roots  of  C(z)  on  the 
unit  circle  and  in  B  two.   Since  roots  on  the  unit  circle  do 
not  contribute  to  the  integral,  in  regions  A  and  C  the 
integral  vanishes.   The  vanishing  of  the  integral  in  region  A 
gives  the  appropriate  range  of  influence  to  the  space  time 
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point  X  =  y  =  t  =  0.  The  line  DD'  shows  the  position  of  the 
possible  discontinuity  plane.  When  a^  =  A  the  section  of 
the  ray  cone  is  given  in  Figure  3.  Now  the  lines  EE'  and  PF' 
niay  also  carry  discontinuities.  The  regions  A,  B  and  C  have 
the  same  general  properties  in  both  figures. 

Parametric  trigonometric  representations  of  the  ray  cone 
have  been  given  elsewhere,^  but  we  may  use  its  description 
as  the  inverse  of  the  pedal  curve  of  the  l/v  locus  to  obtain 
an  algebraic  parametric  representation.  The  locus  of  the  l/v 
curve  is  given  in  rectangluar  coordinates  as: 

(17)     Q(5,n)  =  1-Ca^  ^  t^lne*  n^y  ^l^l^he*  r?)   =  0. 
The  prescription  for  the  ray  cone  then  gives 
(18a)     f  (5  aS||^  +  ^  52|i^)  ,  II 

(18b)     I  {i   SQ^  +  ^  ^^^)   =  H  ■ 

We  may  eliminate     t]     between   (I7)   and   (I8)   and  obtain 

(19a)  0  =  a^A^eS  .  sal^lial  +  A^)!^  +  2a2A2   {^H^ 

M(a„2^A2)'-a2A2)e.(a2  +  A2)f 

=  e((a2A2e'-(a^A^))la2A^)-  ^(Sal^fi   -(a^+A^))        , 
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{19b)  0=  {|)^^a2  a2  e-   (a^  +  A^)) 


-   (a^'*^  a^^e^)V{a^A,^H^  e*a^„^) 


where     i     is  restricted  by 


(20a)  0  <   ^^  <  rain   (i,   i) 

o        o 


or 


2  2 

/      ^^  /I       1  N    ^   i2    ,  ^o   "^  ^o 

(20b)  max(-75-,  -^)   <  ^     <  — ^ — o-  • 

a        A"^     "  a^  A^ 

o       o  o     o 


As  we  shall  not  be  concerned  with  particularly  detailed 
formulae  for  the  Green's  function,  we  shall  not  try  to  get 
the  most  explicit  form  of  an  answer.  But  if  one  were  concerned 
with  an  actual  computation,  rather  than  study  C(z),  a  poly- 
nomial of  fourth  degree  and  with  complex  coefficients,  one  could 
examine  a  real  cubic.   In  the  region  where  the  Green's  function 
is  different  from  zero  we  may  express  the  roots  of  C(z)  in 
terms  of  the  real  parameters  A,  C^   and  r|  as  e  ^  ,   ^  , 

g-i(C+T]j^  and  e~^^~^^.   From  the  standard  relations  between 
roots  and  coefficients  one  can  finally  get  a  real  cubic  with 
roots  cos(20>  cos(i>^+Ti),  and  co3(i>-'n),  and  the  Green's  func- 
tion may  be  written  in  terms  of  these.   We  shall  adopt  a  much 
simpler  expedient;   z.   and  z   are  respectively  the  roots  of 
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C(z)  Inside  and  outside  the  unit  circle,  z,   and  Zp  are 
the  two  on  the  unit  circle. 

We  are  now  almost  ready  to  give  the  explicit  form  of  the 
Green's  function.   But  first  we  shall  make  a  few  comments  on 
the  unphysical  part  of  the  Green's  function,  which  we  shall 

not  include  in  the  formulas.   If  we  define  a  function  h(r',t) 

P    op 
whose  Fourier  transform  is  (cu  -  a  k  )/D,   then  the  unphysical 

^    zt/p^ 
part  of  Ggg,  G^g^   is  given  by  (see  (6d),  (6e)) 


t 


dt'h(r,t' ) 


The  function  h(r,t)   is  just  of  the  type  in  the  physical  Green's 
function,  and  thus  its  time  integral  will  have  the  same  general 
behavior  so  far  as  the  singularity  surfaces  as  the  physical 
Green's  function.   It  is  possible  to  do  the  time  integral 
explicitly,  but  we  shall  not  explore  the  matter  further.   In 
the  expressions  below  we  do  not  include  the  unphysical  parts  of  B. 

We  shall  express  the  solutions  in  terms  of  the  auxilliary 
real  functions: 

r 


TT  r  i 


r+r 


a     -   t^H(t) 

2 24 

'^       2t  r^i  J 


r  2 

Ar.       Z   -H   1 


r^r 


n      -        t^H(t) 
G3 ^"hp 


r+r' 
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where  the  integral  is  over  the  previously  described  contour, 
Thus  we  have: 


^1  = 


2t^H(t)/         "1  ^___,.,___^.^___^ 

7rr'^^'^h^i-^o^^^i-^l^^^i-^2^  "^2o-^i^^^o-^l^^^o-^2^ 


(22a) 


in  A  and  C. 


in  B 


G^+   iG^ 


(22b) 


2t^H(t )  [         ^1    

o 


\ 


0 


V 


in  A  and  C. 


in  B 


/ 


where  the  regions  A,  B,  C  in  (22)  refer  to  Figures  2  and  J>. 

By  reference  to  equations  (6)  and  (15)  we  see  that  three 

fundamental  forms  appear  in  the  physical  solution  corresponding 

to  k  ,  k   and  k  k   in  the  numerator  of  (13).   Corresponding 
■^        y  ^  y 

to  these  forms  we  Introduce  the  Green's  functions  G   ,  G 

XX '   xy 

and  G„^>   so  that 


(23a) 


(23b) 


(23c) 


^xx  =  ^1  +  Gg^^^  ^^  "  G,sin  20 


xy 


G^sin  20   +  G,cos  20 


G   =  Gi  -  G2C0S  20  +   G,sin  20. 
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After  some  manipulation  we  have  finally: 


(2*a)  p  .  (s  ^  -  1^  Sp)((|)  -  A„2)G^  -  ff^  a,^  .  2j  a^p 


-  (*o(\  h  -  \  y-  \  If 


A^s  |-)(Jg  +1  G      ) 


-  ^o  \  ^  ^xy- 


^    o   S  n3/Xx^  «2x^   .  2xr      "2 


(24b)  u  =  (S  I-  -  S   l^)^(J)  -  A^)G   +^G   +^G 
^   ''   X   ^  p  ax    u^  at^  pt^    o-*  XX   72  ^xy   I? 


yy 


*   (%'*o(\  If  -  \   If)-  ^o^p  ly 


a^s   |r)G 
o  u  at ^  xy 


-  (a^  +  A^)S   I-  (J  G   +  1^  G   ) 
^   '^         o^  u^  ay  ^t  xy   t  ''yy^ 


2  .  .2. 
'o 


^''^^      ^  =(-%  i  -  'l^\  h  -   %  M  (^t^'-  ^o)Gxx 


t-   xy   ^^  yy 


^  iF  «xy  -  !^  ^yyh^o^V  i  -^P  lr)(|  ^xy-f  ^yy) 


(2^^)  Bx  =  -\|r(((f)'-o)Gxx-fi^ 


G    +    il)   G 
xy   ^t-*  yy, 


X  y 


-  ^o(\  ly 
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(2i,e)  By  =  -(S   1^  +  S   |^)(((|)  -  a2)G^^ff2  G^^^  (|)  G 


-  ^o'V  h  *  h   IfXf  «xy  *  I  °yy)  '  ^o^.  If  °xy 


Since  the  results  are  so  formidlble  a  very  qualitative 
description  of  the  solution  may  be  useful.  From  (22)  we 
see  that  all  the  functions  in  (24)  are  different  from  zero 
only  in  region  B  of  Figure  1  and  they  become  Infinite  as 
z^   — >  z  ,  or  as  r/t  approaches  the  boundary  of  B,  and 
we  may  consider  the  disturbance  strongest  there.   At  the 
points  D  and  D'   of  Figure  1  the  singularity  in  the  Green's 
function  is  of  a  higher  degree  than  elsewhere,  so  that  loosely 
speaking  the  disturbance  centers  on  D  and  D' .   When 
a  =  A   Figure  2  is  relevant  and  D,D' ,  E,E' ,  F,F'   are 
points  of  highest  singularity. 

If  we  wished  to  perform  the  differentiations  implicit  in 
(24)  we  must  recall  that  the  expressions  are  distributions  and 
derivitives  of  singular  distributions  are  not  Identical  with 
derivitives  of  the  corresponding  functions,  e.g.. 

If  ^  =  (|,)'log  r  =  -^^  -   .6(x)5(y). 
Such  singular  behavior  is  expected  near  D  and  D' . 
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FOOTNOTES 


1.  A  preliminary  account  of  the  results  in  this  report  was 
presented  at  the  December  1959  meeting  of  the  Plasma 
Physics  Division  of  the  American  Physical  Society  at 
Monterey,  Calif. 

2.  See,  e.g.,  the  paper  of  H.  Grad  In  The  Magna to dynamics 

of  Conducting  Fluids,  ed.  D.  Bershader,  Stanford  University 
Press  (iy:?yj,  and  references  cited  there. 

J>.      Distributions  are  used  essentially  as  described  In  An 

Introduction  to  Fourier  Analysis  and  Generalized  Functions 

Dy  M.  J.  Llghthlll,  Cambridge  University  Press  (I958). 

However,  there  are  certain  added  complications  here  not 
present  in  Llghthlll 's  treatment.   In  a  paper  to  appear 
later  we  study  the  fundamental  solution  of  a  general  linear 
first  order  homogeneous  hyperbolic  equation  with  constant 
coefficients  In  two  space  and  one  time  variables.   With 
classical  methods  we  show  the  functions  given  here  are 
derlvltlves  of  ordinary  functions  so  that  the  Integral  of 
the  fundamental  solution  times  a  C"  function  of  compact 
support  is  naturally  defined  by  Integration  by  parts. 

4.  See,  e.g.,  M.  J.  Llghthlll,  op.  clt. 

5.  M.  J.  Llghthlll,  pp  37-38,  op.  clt.  ■ 

6.  A  preliminary  account  of  these  results  was  presented  by 
the  author  at  the  International  Conference  on  Partial 
Differential  Equations  and  ContlnuiAm  Mechanics  at  Madison, 
Wisconsin  June  1959-   A  paper  will  appear  shortly. 

7.  F.  John,  Plane  Waves  and  Spherical  Means,  Chapter  II, 
Intersclence  Press  (1955). 

8.  See  R.  Courant  and  D.  Hllbert,  Methods  of  Mathematical 
Physics,  vol.  2,  Intersclence  Press,  to  appear.   For  the 
general  theory  as  applied  to  these  equations  see  J.  Bazer 
and  0.  Flelschman,  Phys.  Fluids  2,  366  (1959),  and  H.  Grad, 
op.  clt. 

9.  See  J.  Bazer  and  0.  Flelschman,  op.  clt. 
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FIG.    I 

INVERSE   TO  THE  NORMAL 
SPEED   LOCUS 
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FIG.  2 
SECTION  OF  THE    RAY  CONOID  FOR  Qq^Aq 


FIG.  3 
SECTION  OF  THE   RAY  CONOID  FOR   ao=Ao 
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